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We analyze recent heat capacity measnrements in herbertsmithite ZnCu 3 (OH) 6 Cl 2 single crystal 
samples subjected to strong magnetic fields. We show that the temperature dependence of specific 
heat Cmag formed by quantum spin liquid at different magnetic fields B resembles the electronic 
heat capacity Cei of the HF metal YbRh 2 Si 2 . We demonstrate that the spinon effective mass 
^mag ^ Cmag fT exhibits a Scaling behavior like that of CeijT. We also show that the recent 
measurements of Cmag are compatible with those obtained on powder samples. These observations 
allow us to conclude that ZnCu 3 (OH) 6 Cl 2 holds a stable strongly correlated quantum spin liquid, 
and a possible gap in the spectra of spinon excitations is absent even under the application of very 
high magnetic fields of 18 T. 

PACS numbers: 75.40.Gb, 71.27.+a, 71.10.Hf 


Quantum spin liquids (QSLs) are promising new 
phases, where exotic quantum states of matter could be 
realized. Although much theoretical effort has been de¬ 
voted to understand their physical nature, the question is 
still far from its complete clarification. Generally speak¬ 
ing, the QSL is a quantum state, formed with hypothetic 
particles like fermionic spinons carrying spin 1/2 and no 
charge. A number of QSLs with various types of ground 
states are proposed |ll4l8l| but the lack of real materials 
possessing them obscure the underlying physical mecha¬ 
nism. On the other hand, one needs a real theory that 
plays important role in the understanding and interpret¬ 
ing accessible experimental facts [igl - E^ . Measurements 
on magnetic insulators with geometrical frustration pro¬ 
duce important experimental data shedding light on the 
nature of spinon composing QSL. Recent measurements 
indicate that the insulator ZnCu 3 (OH) 6 Cl 2 (herbert¬ 
smithite) is very likely to be the first promising candidate 
to host a QSL in real bulk materials. In herbertsmithite, 
the dynamic magnetic susceptibility shows that at low 
temperatures quasiparticle excitations, or spinons, form 
a continuum, and populate an approximately flat band 
crossing the Fermi level [Sim- At the same time, our 
analysis of herbertsmithite thermodynamic properties al¬ 
lows us to reveal their scaling behavior. The above results 
demonstrate that the properties of the herbertsmithite 
are similar to those of heavy-fermion (HF) metals. Thus, 
it can be viewed as a new type of strongly correlated HF 
electrical insulator exhibiting p ro p erties of HF metals 
but resisting electric current [20|-l22|. 


The development of this concept, however, faces funda¬ 
mental problems still remaining to be resolved. The first 
one is that the experimental data are taken in measure¬ 
ments on herbertsmithite powder samples 
As a result, both out-of-plane magnetic defects and site 


disorder between the Gu and Zn ions can strongly change 
the behavior of the bulk susceptibility or even alter its 
low-t emp erature scaling at weak magnetic fields, see e.g. 
Ref. [1^. The second problem is based on experimental 
data showing that QSL in the herbertsmithite is unstable 
against the application of external magnetic fields. This 
instability is represented by the emergence of so-called 
spin-solid phase, which, in turn, is separated from the 
spin-liquid one by a gap induced by an applied magnetic 
field [23|. Thus, the experimental observations of QSL 
remain scattered and need regimentation. In our opin¬ 
ion, this can be partially achieved by presenting a reliable 
interpretation of the recent measurements of specific heat 
C in high magnetic fields performed on herbertsmithite 
single crystal samples 2J, [2^ . 


In this Communication we analyze the recent heat ca¬ 
pacity C measurements in strong magnetic fields on her¬ 
bertsmithite. We interpret the above measurements in 
terms of QSL spinon contribution. To obtain the spinon 
- determined partial specific heat Cmag, we subtract the 
phonon and Schottky effect contributions from the to¬ 
tal heat capacity C. It turns out that obtained Cmag 
as a function of temperature T at fixed magnetic field 
B behaves very similarly to the electronic specific heat 
Cel of the HF metal YbRh 2 Si 2 . This observation allows 
us to conclude that spinons form the Fermi sphere with 
the Fermi momentum pp, while the gap in the spectra 
of spinon excitations and the spin-solid phase are ab¬ 
sent. We demonstrate that the effective mass of spinon 
^mag ^ Cmag IT exhibits a scaling behavior resembling 
that of CeijT■ We also show that new measurements of 
Cmag On single crystal samples are compatible with those 
obtained on powder samples. 

The S' = 1/2 spins of the Cu+^ ions occupy the posi¬ 
tions in a highly symmetric kagome lattice. As spins on 
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FIG. 1: (color online). The heat capacity measurements on 
single crystal of ZnCu^OH) 6 Cl 2 at B = 0 (squares) and B — 
14 T (triangles) [13, [Is] • The thick blue line corresponds to 
our theoretical approximation m with fitting parameters a 
and b. The lines represent C = and C = bT^, see 

legend. 


the above lattice reside in a highly symmetric structure 
of corner-sharing triangles, the ground state energy does 
not depend on the spins configuration. As a result, spins 
located at the kagome hexagon, composed of the two tri¬ 
angles, form a (frustrated) pattern that is even more frus¬ 
trated than the triangular lattice considered by Anderson 
[i^ . As a result, the kagome lattice has a flat topologi¬ 
cally protected branch of the spectrum with zero excita¬ 
tion energy Therefore, the fermion condensation 


quantum phase transition (FCQPT) can be considered as 
a quantum critical point (QCP) of the ZnCu 3 (OH) 6 Cl 2 
QSL In that case, we expect that at low temper¬ 

atures the heat capacity C of herbertsmithite becomes a 
function of the applied magnetic field, as it is seen from 
Fig. □ We propose that mentioned QSL is composed 
of chargeless fermions called spinons with the effective 
mass Latter fermions with spin ct = 1/2 occupy 

the corresponding Fermi sphere with the Fermi momen¬ 
tum pp, and form the excitation spectrum typical for HF 
liquid located near FCQPT, while spinons represent HF 
quasiparticles of deconfined QSL. Thus, the ground state 
energy E[n) is given by the Landau functional depend¬ 
ing on the spinon distribution function no-(p), where p is 
the momentum. Near FCQPT point, the effective mass 
Ml 


is governed by the Landau equation [291 - 1311 
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Here we have rewritten the spinon distribution function 
as i5ncr(p) = n^rip, T, B)—ncr(p, T = 0, H = 0). The Lan¬ 
dau interaction F is defined by the fact that the system 


has to be at FCQPT. Thus, this interaction brings the 
system to FCQPT point, where Fermi surface alters its 
topology so that the effective mass acquires temperature 
and field dependence At this point, the term 

l/Mll^g{T = 0, B = 0) vanishes and Eq. ([T|) becomes ho¬ 
mogeneous. It can then be solved analytically. At H = 0, 
the effective mass depends on T demonstrating the NFL 
(non-Fermi liquid) behavior 

M*^^giT) ^ aTT-^/\ ( 2 ) 

At finite T, the application of magnetic field B drives the 
system to Landau Fermi liquid (LFL) region with 

M*^,g{B) c. aBH-2/3. (3) 


At finite B and T near FCQPT, the solutions of Eq. 
o Ml^ag{B,T) can be well approximated by a sim¬ 
ple universal interpolating function. The interpolation 
occurs between the LEL (M*(T) oc const) and NEL 
{Ml,,g{T) oc r-2/3) regions. As it is seen from Eq. 

under the application of magnetic field M^^^g be¬ 
comes finite, and at low temperatures the system demon¬ 
strates the LFL behavior Cmag{T)/T oc Ml^^g{T) ~ 
Fllnagi'P' = 0) -|- oiT^. Then, as it is seen from Fig. 
[21 at increasing temperatures M^^^g oc CmagjT grows, 
and enters a transition regime, reaching its maximum 
OC {Cmag{T)/T) max 3jt T — -^maxi ^Vltrh SubsCCJUGHt 
diminishing given by Eq. 

To reveal a scaling behavior, we introduce the nor¬ 
malized effective mass and the normalized temper¬ 
ature Tn dividing the effective mass M^^^g by its maxi¬ 
mal values, and temperature T by Tmax at which 

the maximum occurs. The normalized effective mass 
= M* as a function of the normalized tem¬ 

perature y = Tp[ = T/Tmax is given by the interpolat¬ 
ing function that approximates the solution of Eq. o 
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Here cq = (I-I-C2)/(1 + Ci), c\ and C 2 are fitting parame¬ 
ters, Cl and C 2 are fitting parameters, parameterizing the 
Landau interaction. Magnetic field B enters Eq. ([T]) only 
in the combination pbB/ ksT, making ksTmax — ysB 
where kp is Boltzmann constant and /is is Bohr magne¬ 
ton 3^, 321. Thus, Eq. (j?]) exhibits the scaling behavior, 
which is the intrinsic property of Eq. (|T|) at FCQPT: M'^ 
at different magnetic fields merges into a single curve as 
a function of the variable y. In what follows we use Eq. 
(|4|) to clarify our calculations based on Eq. (HD. 

To analyze the spinon specific heat Cmag{T, B), we 
have to exclude the contribution coming form lattice 
(phonons) that is independent of H. To separate Cmag 
contribution, we approximate C{T,B = 0) by the func¬ 
tion 


C(T) = ar^/3 


( 5 ) 
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FIG. 2: (color online). Panel a: The specific heat CmagjT 
of ZnCuafOHleCb measured on powder and single crys¬ 
tal samples [Ij, [^. (Cmag/T) is displayed as a fnnction of 
temperature T for fields B shown in the legends. Panel 
b: The specific heat Csi/T extracted from measurements on 
the HF metal YbRh 2 Si 2 [Ij exhibits the same behavior as 
Cmag/T. The illustrative values of the maxima {Cmag/T)raax 
and (Cei/T) max of both CmagiT and Ce.ilT and the corre¬ 
sponding Tmax are also shown in panels a and b. 


Here the first term proportional to M^^g{T) oc 
presents the QSL contribution when it exhibits the NFL 
behavior, as it follows from Eq. while the second 
one is determined by the lattice (phonon) contribution. 
It is seen from Fig. [T] that the approximation dSD is 
valid in a wide temperature range. At low temperatures 
T < 4 K the heat capacity is mainly formed by spinons, 
while the phonon contribution prevails at higher temper¬ 
atures, making its main contribution at T > 10 K. It is 
seen, that at diminishing temperatures the LFL behav¬ 
ior, C oc T, starts to form. We note, that the Schot- 
tky anomaly was removed, as it is usually done, see e.g. 
[2J, [ 23 . Upon comparing the data on the heat capac¬ 
ity taken at H = 9 and H = 14 T shown in Fig. [1] 



FIG. 3: (color online). Schematic T — B phase diagram of SC- 
QSL with magr^etic field as control parameter. The vertical 
and horizontal arrows show LFL-NFL and reverse transitions 
at fixed B and T, respectively. The hatched areas in the 
inset and main panel shown by the arrows represent the tran¬ 
sition region at Tmax(H). The solid line in the hatched area 
represents the function rmax(H) ~ Tm{B). The functions 
Tm{B) shown by two-headed arrows define the total widths 
of the NFL state and the transition area. The inset shows a 
plot of the normalized effective mass versus the normal¬ 
ized temperature Tn- Transition region, where M* reaches 
its maximum at Tn = T /Tmax = 1, is shown by the 

arrow. 


one concludes that C strongly depends on the magnetic 
field. However, at raising temperatures T > 10 K the B- 
dependence vanishes, for the contribution coming from 
phonons prevails, then, the system enters the transition 
regime at ksT ~ TbB, and the influence of magnetic 
field becomes weak, as shown in the T — B phase dia¬ 
gram, Fig. m 

The obtained heat capacity of QSL, Cmag/T = {C — 
bT^)/T, is displayed in Fig. El panel a. It follows from 
panel a, that the results obtained on different kinds 
(powder and single crystal) of samples and in different 


iiwie c 

Si 


measurements [J, [^, 1241 l25l| exhibit similar properties. 
To show that the behavior of Cmag/T reported in panel 
a is of generic character, we note that Cmag/T oc 
behaves like Cei/T shown in panel b. Indeed, as seen 
from Fig. El the maximum structure in Cmag/T and 
Cei/T at temperature Tmax appears under the applica¬ 
tion of magnetic field B and Tmax shifts to higher T 
as B is increased. The value of the maximum struc¬ 
ture is saturated towards lower temperatures decreasing 
at elevated magnetic fields. It is worthy to note, that 
Cmag/T ^ Cei/T, as it is Seen from the panels a and 
b. Thus, QSL of ZnCu 3 (OH) 6 Cl 2 can be viewed as a 
strongly correlated quantum spin liquid (SCQSL) that 
behaves like that of the HF metal YbRh 2 Si 2 , exhibiting 
a gapless behavior in high magnetic fields up to 18 T. 

We now construct the schematic phase diagram of 
ZnCu 3 (OH) 6 Cl 2 , reported in Fig. El At T = 0 and 
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FIG. 4: (color online). Panel a: The normalized specific 
heat {Cmag/T)N = of ZnCu 3 (OH) 6 Cl 2 versus normal¬ 
ized temperature Tjv as a function of B shown in the legends. 
(Cmag /T)n is extracted from the data [i, m, m shown in 
Fig- m panel a. Solid curve represents our calculations, and 
traces the scaling behavior of the effective mass. Panel b: The 
normalized specific heat {Cei/T)^ = of YbRh 2 Si 2 versus 
normalized temperature Tjv as a function of B at low (upper 
box symbols) and high (lower box symbols, corresponding 
to those shown in Fig. [2] panel b) magnetic fields, extracted 
from the specific heat measurements on the YbRh 2 Si 2 [33.[35l|. 
The low-field calculations are depicted by the dotted line trac¬ 
ing the scaling behavior of M)). The high-field calculations 
(solid line) are performed for magnetic fields when the quasi¬ 
particle band becomes fully polarized [3^ . 

B = 0 the system is approximately located at QCP of 
FCQPT without tuning. Both magnetic field B and tem¬ 
perature T play the role of the control parameters, shift¬ 
ing the system from its QCP and driving it from the NFL 
to LFL regions as shown by the vertical and horizontal 
arrows. At fixed temperatures the increase of B leads 
the system along the horizontal arrow from NFL to LFL 
region. This behavior is seen from Fig. [2J panels a and b: 



FIG. 5: (color online). Panel a: The maxima (Cmag/r)max 
of Cmag It versus magnetic field B are shown by the solid 
squares, see Fig. [2] panel a. The solid curve is approximated 
by Mm^^(B) = dB~^^^, see Eq. ([3ll, where di is a fitting 
parameter. Panel b: The temperatures rmax(B), at which 
the maxima of (Cmag/T) are located, see Fig. [2 panel a. 
The solid straight line represents the function Tmax = d 2 B, d 2 
being a parameter. Td shown by the arrow is the temperature 
at which the straight line Tnax starts to deviate from the 
experimental data, as shown by the dashed line. 


The application of B shifts Tmax to higher temperatures. 
On the contrary, at fixed magnetic field and increasing 
temperatures the system moves along the vertical arrow 
from the LFL to NFL region. The inset to Fig. |3] demon¬ 
strates the behavior of the normalized effective mass M^ 
versus the normalized temperature Tjv that follows from 
Eq. (|4]). It is seen that the temperature range Tn ~ 1 
represents the transition region between the LFL regime 
with almost constant effective mass and the NFL behav¬ 
ior, having the T“2/3 dependence. It follows from Eq. o 
and Fig. [3]that the width of the transition region is given 
by Tm C!C T^a^[B) (xT (X B. The phase diagram in Fig. 
[3] demonstrates that the properties of ZnCu.- 3 (O H)fi CE 
are close to those of the HF metal YbRh 2 Si 2 [Sd. l3l|. 

Figure m panel a, reports the normalized specific heat 
{Cmag/T)N = versus normalized temperature 
as a function of B. It is seen that at low T/v < 0.1, 
{Cmag/T)N — 0.4. This value is determined by the QSL 
polarization in magnetic fields, and coincides with that of 
{Cei/T)N = obtained on YbRh 2 Si 2 . In both panels a 
and b, our calculations are shown by the same solid curve. 
Note that at low temperatures and magnetic fields, when 
the polarization is negligible, {Cei/T) at ~ 0.9, as it is seen 
from panel b of Fig. [H 

In Fig. 0 panel a, the solid squares denote the val¬ 
ues of the maxima [Cmag/T)^a_^[B), taken from Fig. O 
panel a. In panel b, the corresponding values of 'rmax(.S) 
are displayed versus magnetic field B. It is seen that the 
agreement between the theory and experiment is good. 
In panel b, the solid straight line represents the function 
Tmax(-B) c)c c?2 -B- R is Seen that at relatively high tem¬ 
peratures the experimental data are well approximated 
by the straight line, while at T < there is deviation. 
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This fact evidences that SCQSL of herbertsmithite is not 
exactly placed at FCQPT, for in that case the approxi¬ 
mation Tinax oc d 2 B would be good at low temperatures. 
This conclusion coincides with general properties of the 
phase diagrams of HF metals 37| and with experimental 
data, showing that at low temperatures C demonstrates 
the LFL behavior [ 1 , 24, 2^, see Fig. [T]as well. 

In summary, based on the recent experimental facts, 
we have shown that ZnCu 3 (OH) 6 Cl 2 can be viewed as a 
strongly correlated Fermi system whose thermodynamics 
is defined by SCQSL located at FCQPT. Our calcula¬ 
tions of the heat capacity Cmag/T are in good agree¬ 
ment with the experimental data and the scaling be¬ 
havior of Cmag/T coincides with that observed in the 
HF metal YbRh 2 Si 2 up to high magnetic fields of 18 T. 
Thus, spinons forming SCQSL remain itinerant even in 
high magnetic fields. We also have shown that recent 
measurements of Cmag on single crystal samples are con¬ 
sistent with those obtained on powder samples. 
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